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ABSTRACT 



In most shipboard CIC's true motion plots of own ship 
and/or surface/subsurface contact motions are performed 
with the aid of electromechanical devices. These are 
inaccurate and have a variety of drawbacks. 

This thesis describes a fast, labor-saving, bearing 
only radar tracking solution utilizing various microprocessors. 
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I. 



INTRODUCTION 



Passive localization and tracking techniques are of 
interest in a variety of microcomputer based surface/ 
subsurface radar applications. 

Because a significant portion of the cost of building 
or maintaining a warship is the electronics in its sensor 
and weapon system, and because a great amount of time and 
manpower is employed in performing simple but important 
tasks, microcomputers offer the potential to: 

(1) Reduce the cost of digital systems, 

(2) Perform some complex functions at remote stations, 
relieving the congestion at larger central computer 
facilities, and 

(3) Perform functions currently handled by watch personnel, 
thus reducing the manning requirements of watch 
sections . 

An example of this is the problem of manual tracking of radar 
contacts and the solution of Maneuvering Board Problems. 

It is the purpose of this study to demonstrate than an 
alternate approach to the solution of the problems mentioned 
above can be developed and implemented by using a micro- 
computer system, while at the same time maintaining a human 
engineered user interface. 
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II. FORMULATION OF THE PROBLEM 



A. PREFACE 

For Naval ships not equipped with Naval Tactical Data 
Systems (NTDS) , operations performed by the Combat Infor- 
mation Center (CIC) during a normal peacetime watch include 
manual tracking of radar contacts and solution of maneuvering 
board problems. The environment emulated by the Intel 
Intellec Microcomputer Development System (MDS) provides a 
data subset of the NTDS. The subset of information is 
provided via radar interface to the MDS System. 

The interface between the microcomputer and the surface 
radar was not constructed to due cost, time and manpower 
constraints . 

The scenario of the study was limited to the development 
of a Kalman Filter bearing-only tracking algorithm for the 
Intel Intellec MDS microcomputer system and Intel SBC-310 
High-Speed Mathematics Unit Performance. Simulation of the 
circuitry is necessary to interface with the AN/SPS-10 
surface search radar. 

The parameters for the AN/SPS-10 surface search radar 
were used for computations requiring specific radar 
parameters. This radar was chosen because all information 
regarding the radar is available at the Naval Postgraduate 
School . 
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A Kalman filter was selected because it minimizes the 



estimation error in a well defined statistical sense and 
most complex calculations of the tracking filters typically 
used for bearing-only radar tracking that is easily 
computable with a microcomputer system. 

The system input/output parameters are shown in 
Figure 1. 

B. TYPICAL COMBAT INFORMATION CENTER OPERATIONS 

During normal peacetime steaming, the CIC watch team 
may consist of from two to ten or even more personnel, 
depending on the size of the ship as well as on the complexity 
of the equipment being used. 

Among the problems that are normally solved by CIC 
personnel^ special mention needs to be made of those of 
plotting contacts and the determination of parameters such as 
course, speed and closest point of approach (CPA) of those 
contacts. This is a tedious and error-prone task; it often 
requires most of the time and effort of the CIC team and it 
is vitally important to the safety of the ship. This had 
led to the installation of equipment to reduce the amount 
of workload in the CIC while at the same time improving the 
reliability of the constant information provided to the 
bridge. This equipment includes dead reckoning devices and 
NC2 plotter. 

1 . Maneuvering Board Plotting Sheets 

The maneuvering Board Plotting Sheets (H.O. 2665-10) 
have been prepared in order to facilitate the solution of 
a ship's relative movement problem. 



Oontact 

Bearing 


SYSTEM 


Expected Contact Course 


AN/SPS-10 

radar Oontact 

Range 


Expected Oontact Speed 


FROM GYRO Ship 

COMPASS Course 


Expected Oontact Range 


AND 

PITOMETER Own Ship 

ICG Speed 


Expected Contact Bearing 






Time 









INPUTS 



OUTPUTS 



Figure 1. SYSTEM INPUT/OUTPUT PARAMETERS 
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The primary responsibility of CIC is to provide 
information and recommendations on the tactical situations, 
such as contact course, speed, CPA information, are defined 
using the "Maneuvering Board" plotting sheets. But the use 
of the Maneuvering Board requires some pracice, also it 
requires the complete attention of one person during the 
CIC operation. 

2 . Dead-Reckoning Equipment 

This equipment maintains a continuous, up to the 
minute, geographic plot of own ship's plot in the CIC. 

The Dead-Reckoning System consists of the following 
basic components: (1) Dead-Reckoning Analyzer (DRA) ; 

(2) Dead-Reckoning Indicator (DRI) ; and (3) Dead-Reckoning 
Time (DRT) . Course and speed inputs of own ship are fed 
into the DRA from Gyrocompass and pitometer-log and then to 
the DRT, where they cause a movable source of light to trace 
the ship-track continuously. 

In all systems as mentioned above the target's rela- 
tive position measurements are obtained from radar repeater 
and plotted on the systems which are read to be used. The 
range measurement accuracy is at best ± 20 yards and 
bearing measurement accuracy is ± 1 degrees with a well 
calibrated repeater and a well trained operator. The own 
ship speed and course are used' to determine the own ship's 
velocity vector. The course and speed of own ship may vary 
± 1 degree and ± 1 knot depending on the helmsman, weather 
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conditions and the ship instruction. The combined result 
and the manual solution will cause ± 5 degrees in course 
and ± 3 knots in speed of the target information. Also 
the major disadvantage of the manual solution, in addition 
to its poor accuracy, is the time required to obtain the 
solution and the inability to obtain solutions of either 
the target or own ship maneuvers between radar position 
measurements, typically 3 or more minutes apart. 

C . MICROPROCESSOR TECHNOLOGY 

Recent advances in large scale integration (LSI) semi- 
conductor process technology have made possible substantial 
reductions in the cost and size of digital logic circuits. 
Microprocessors represent a very remarkable achievement 
of engineering ingenuity and industrial know-how at their 
best. 

Since 1973, the year in which Intel Corporation shipped 
the first 8080, 8-bit N-channel microprocessor, a number of 
manufacturers have developed similar products and because 
of this competition prices have been drastically lowered. 
With hardware and software system development costs as they 
are, much micro work is done on a custom basis. The primary 
advantage of microcomputer based systems is to have 
specialized and dedicated equipment solving specific 
problems . 

The microprocessor families are: (1) 4-bit machine; 

(2) 8-bit machine; (3) 16-bit machine; and (4) bit-slices 
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architecture. 16-bit microprocessor is the best solution 
for our concern of tracking problem. But Intel's 8080 was 
selected as an 8-bit microprocessor and programmed for 
double precision arithmetic to perform Kalman Filter 
algorithms, because of its lower cost and availability at 
the Naval Postgraduate School with all its hardware and 
software support devices. 

D . SCENARIO 

In order for a ship to maintain operational readiness 
as a unit of a task force, it must be aware of the current 
operational environment. The operational environment is 
defined as the surface/subsurface contact profiles in the 
geographic area of interest. The contact profile consists 
of friendly, hostile, and unknown contacts with associative 
contact characteristics. Contact characteristics are such 
measurements from radar, range and bearing to accurately 
determine the course and speed of a radar contact. If the 
radar measurements were perfectly accurate and the target 
had zero acceleration, then the tracking problem would be 
velocity vector problem as shown in Figure 2, velocity 
vector and contact relative motion plot as shown in Figure 3. 
But the radar bearing and range measurements on any scan have 
errors which are a function of the radar system and the 
received signal to noise ratio. The measurement errors will 
be discussed after defining the Kalman Filter parameters. 
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Figure 2. GEOGRAPHIC PLOT OF CONTACT AND OWN SHIP MOTION 
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Figure 3. VELOCITY VECTOR AND CONTACT RELATIVE MOTION PLOT 
(Own Ship assumed at the origin.) 
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The contact motion on a X,Y coordinate system can be 
approximated by linear dynamic system. The target dynamics 
can be represented in a matrix form as shown below to define 
the sample data system such as a search radar where the 
measurements are taken once each antenna rotation 
(Reference 13) 
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vehicle state vector at scan n, 

E/W position of the vehicle, 

E/W velocity of the vehicle, 

N/S position of the vehicle, 

N/S velocity of the vehicle, 

interval between observations, 

acceleration acting upon the vehicle from 
scan (n) to (n+1) . 
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A Kalman filter is a common optimal filtering technique 
for estimating the state of a linear system. It is selected 
to implement the tracking problem because it is computa- 
tionally the most demanding technique in our concern. The 
details of the Kalman filter will be analyzed in Chapter IV. 
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III. SYSTEM DESIGN CONSIDERATIONS 



A. HARDWARE DESIGN CONSIDERATIONS 

Microprocessors, in general, are not complete computers, 
but Central Processor Units (CPU) implemented with one to 
ten' large-scale integrated-circuit chips. Large-scale 
integration (LSI) chips are comprised of 1,000 or more 
gates; many LSI chips hold over 6,000 gates or a "complete 
central processor." 

The word "micro-computer" refers to a small stored 
program computer comprising memory and input/output circuits 
together with a microprocessor CPU. 

Interface design and programming make the mass-produced 
microcomputers into new special-purpose dedicated machines. 

The Intellec Microcomputer Development System (MDS) is 
a complete, coordinated computer system designed around 
Intel's 8080 microprocessor. The MDS has a 2 microsecond 
instruction cycle, a repertoire of 72 powerful instructions, 
unlimited subroutine nesting, and a versatile interrupt 
scheme . 

The MDS was selected because of its 8080 microprocessor, 
its memory capabilities, its interrupt mechanism, and its 
extended I/O capabilities. Also, the selection of the 
Intel Microcomputer Development System (MDS) was made for 
this theis because of its immediate availability at the 
Naval Postgraduate School. 
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A Datamedia Elite 2500 Video Terminal was chosen to 



present alphanumeric information because its features 
include; (1) Editing and roll operation modes; (2) 50 to 
3600 baud programmable speed transmission; (3) protected 
fields, (4) computer derived or high light field (blink) ; 

(5) addressable cursor. 

Because of the computation of Kalman filter algorithms, 
floating-point calculation required, an SBC-310 high-speed 
math unit, developed by Intel Corporation, was incorporated. 

In performing high speed mathematical functions, the Math 
Unit acts as an intelligent processor, performing a reper- 
toire of 14 arithmetic functions and at least an Oder of 
magnitude faster than comparable software routines. Why 
it is to selected to perform the algorithms will be discussed 
in Section C of this Chapter, and Figure 4 shows the final 
configuration of the equipment used in developing this 
thesis. Also, start-up procedures for the system were listed 
in Appendix B. 

B. SOFTWARE DESIGN CONSIDERATIONS 

1 . Language Selection 

Choice of the computer language to be used was 
partially dependent upon the hardware configuration. 

PL/M-80 was preferred as a programming language to be used, 
after selecting MDS hardware system, to the assembly language 
because of its ease as programming and immediate availability. 

PL/M-80 is a high level language developed by Intel 
Corporation and esigned especially for system and applications 
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Figure 4. EQUIPMENT CONFIGURATION 

prograrnming for the Intel 8080 microprocessor. It is a 
block structxired language which speeds the microprocessor 
program development by relieving the programmer of the 
tedious task of memory and register management. 

ISIS-II disk operating system also developed by Intel 
Corporation for the Intellec MDS system has a resident 
PL/M-80 compiler which has been very useful during the 
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implementation, debugging and testing the program. The 
PL/M-80 compiler accepts the statements as input and produces 
a machine-code program module as output. 

PL/M statements are divided into two basic categories 

1 . DECLARE and PROCEDURE statements . 

2. EXECUTABLE statements, which are all other than declare 
and procedure statements. 

EXECUTABLE statements cause machine code to be generated. 
DECLARE statements cause variables to be defined associated 
with objects and PROCEDURE capability permits program to 
be constructed in a modular fashion, which has niamerous 
advantages including efficiency of coding, readability of 
programs, ease of debugging and possibility of using the 
same procedure in more than one program (Ref. 6)'. 

2 . Floating Point Arithmetic 

Floating point arithmetic is required because range 
of numbers to be represented is large and sometimes unpre- 
dictable, especially in the Kalman filter algorithm. In 
this algorithm when the time concern to use a floating point 
arithmetic is disadvantageous, because it needs additional 
subroutines or hardware, but when the accuracy concern we 
must use it in our programming. 

In selecting floating-point formats, the primary 
considerations are word size and the radix of the arithmetic. 
In much of the literature, 24-bit word size is regarded as 
too small for scientific computation. The Intel standard 
for floating point arithmetic has been adopted to apply to 
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all general purpose products including software systems, 
and components. Such as FPAL (Floating-Point arithmetic 
library) for the 8080 and Math Board (SBC-310) made of series 
3000 bit slices have already been produced using this 
standard. 

As the floating-point arithmetic formats concern 
there are two standard formats which have been accepted 
from Intel (Ref. 12); (a) 32-bit short precision word; 

(b) 62-bit long precision word. The 32-bit short precision 
word format with binary radix was chosen because of the hard- 
ware consideration by selectin the MDS system and 8080 
microprocessor with SBC-310 high speed mathematics unit was 
the fastest and available at the Naval Postgraduate School. 

Comparison will be done next with its compatibles: 

(1) AM9511 Arithmetic processing unit (Advanced 
Micro Devices, Inc.) 

(2) 8080/8085 Floating-Point Arithmetic Library (FPAL) 
(Intel Corporation) 

(3) 24-bit special floating point (Ref. 11) 

(4) SBC 310 High Speed Mathematic Unit (Intel 
Corporation) . 

a. AM9511 Arithmetic Processing Unit 

The AM9511 is a monolithic MOS/LSI device which 
provides high speed, 16 and 32 bit fixed and 32 bit floating- 
point arithmetic, plus a group of transcendental derived 
functions, control and conversion commands. 

The AM9511 APU is contained within a single 
dual-in-line 24 pin chip as described in detail in Ref. 3, 
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The distinctive characteristics of the AM9511 



APU include; 

- Fixed-point 16 and 32 bit operation 
Floating-point 32 bit operation 

- Binary data formats 

Basic add, subtract, multiply and divide 
Trigonometric sine, cosine, tangent 

Inverse trigonometric arcsine, arccosine, arctangent 
Square roots 

Logarithms base on 10 and e 
Exponentiation 

Float-to-f ixed and f ixed-to-float conversions 
stack oriented operand storage 
DMA or programmed I/O data transfers 
General purpose 8-bit data interface. 

Data Formats ; 

All fixed-point operands and results are 
represented as binary two's complement integer values. The 
16 bit format can express nximbers with a range of -32,768 
to +32,767. The 32 bit format can express numbers with a 
range of -2,147,483,648 to +2,147,483,647. 

The floating point format uses a 32 bit word 
width as shown in Figure 5. The most significant bit (bit 31) 
indicates the sign of mantissa. The next seven bits form 
the exponent and remaining 24 bits form the mantissa value. 
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EXPONENT SIGN 

EXPONENT MANTISSA 

jslsl I 1 i I 

3130 24 16 18 0 



Figure 5. FLOATING-POINT FORMAT FOR AM9511 APU 

The exponent of the base 2 is an unbiased two's 
complement number with a range of -64 to +63. The mantissa 
is a sign-magnitude number with an assumed binary point 
just to the left of the most significant mantissa bit (bit 2) . 
All floating point values must be normalized which m.akes 
bit 23 always equal to 1 except when representing a value of 
zero. The number zero is represented with binary zeros in 
all 32 bit positions. 

The operation of the APU is described in Ref. 3 
in detail with the command formats and their descriptions. 

The AM9511 interested commands, its clock cycles and 
typical execution times for three different clock frequencies 
such as 2 MHz, 3 MHz and 4 MHz are shown in Table I. 
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TYPICAL EXECL7TION TIMES 


OCMIAND 

MNEMONIC 


CLOCK CYCLES 


AM9511 
(2 MHz) 


AM9511-1 
(3 MHz) 


AM9511-4 
(4 MHz) 


32-BIT FIXED POINT OPERATIONS 


ADD 

(DADD) 


22-22 


10 


7 


5 


SUBTRACT 

(DSUB) 


38-40 


19 


13 


10 


MULTIPLY, LCI-ZER 
(DMUL) 


184-210 


97 


65 


49 


MULTIPLY, UPPER 
(DMUU) 


182-218 


91 


61 


46 


DIVIDE 

(DDIV) 


186-210 


98 


66 


50 


32-BIT FLOATING-POINT PRIMARY OPERATIONS 


ADD 

(FADD) 


54-368 


27 


18 


14 


SUBTRACT 

(FSUB) 


70-370 


35 


24 


18 


MULTIPLY 

(EMUL) 


146-168 


73 


49 


37 


DIVIDE 

(FDIV) 


154-184 


77 


52 


38 


32-bit FLOATING-POINT DEPTVED OPERATIONS 


SQUARERDOT 

(SORT) 


782-870 


391 


262 


200 


LOG BASE 10 
(LOG) 


4474-7132 


2237 


1493 


1125 


LOG BASE e 
(IW) 


4298-6956 


2149 


1435 


1098 



Note: All time values are specified in microseconds. 

Table I. AM9511 ARITHMETIC PROCESSING UNIT COMMAND 
EXECUTION TIMES 
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TYPICAL EXECUTIOJ TIMES 


COMMAND 

MNEMCJIIC 


CLOCK CYCLES 


AM9511 
(2 MHz) 


AM9511-1 
(3 MHz) 


AM9511-4 
(4 MHz) 


CONTINUATION OF 32-BIT FLQATING-POINr DERIVED OPERATIONS 


SINE 

(SIN) 


3786-4808 


1898 


1266 


950 


COSINE 

(COS) 


3840-4878 


1920 


1280 


960 


TANGENT 

(TAN) 


4894-5886 


2447 


1635 


1225 


ARCSINE 

(ASIN) 


6230-7938 


3115 


2080 


1560 


ARCCOSINE 

(AOOS) 


6304-8284 


3152 


2105 


1580 


ARCTANGENT 

(ATAN) 


4992-6536 


2496 


1665 


1250 


EXPCNENT 

(EXP) 


3784-4878 


1887 


1265 


949 


POWER 

(PWR) 


8280-12032 


4145 


2765 


2075 


32-BIT DATA AND STACK MANL 


PULATICN OPERATIONS 


FIXED-TO-FLQAT 

CONVERSION 

(FLTD) 


56-342 


28 


19 


14 


FIOAT-TO-FIXED 

CONVERSION 

(FXTD) 


20-336 


45 


30 


23 


float point sign 
1 CHANGE 
(CHSF) 


16-20 


8 ! 


6 


4 


EXCHA^JGE 

(XCHF) 


26 


13 


9 


7 



Note: All time values are specified in microseconds. 



Table 1 (Continuation) . AM9511 ARITHMETIC PROCESSING 

UNIT COMMAND EXECUTION TIMES 
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b. 8080/8085 Floating-Point Arithmetic Library (FPAL) 
The FPAL procedures reside in object code module 
form in the library "FPAL. LIB" on the Intellec ISIS-II 
System diskette. They are self contained and can be used 
in component, OEM-board, or Intellec Microcomputer, Develop- 
ment System. It contains basic floating point subroutines 
and functions (referred to generally as "procedures") . 

The operations provided are: 

- Addition, subtraction, multiplication, division 

- Value comparison 

- Conversion between decimal and binary floating point 
and 32 bit signed integer formats 

- A default error handler subroutine. 

For all operations single precision format is used and in 
addition to these operations, a number of procedures are 
provided to deal with the Floating Point Record (FPR) . This 
is a reserved, 18 byte work area used to collect status and 
error information, and as an accumulator for intermediate 
results. The procedures supporting the FPR perform FPR 
initialization, change error recovery options, check the 
contents of FPR fields, and pass numbers between the FPR 
and memory. 

The FPAL also includes a default error-handler 

t 

subroutine. This subroutine is called when an invalid number 
is used in a floating-point operation or if overflow, 
underflow, or division by zero are not handled by an 
arithmetic subroutine. We may also write our own error 
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handler, so long as it conforms to the formats described 
in Reference 2. 

The FPAL can be used by assembly langugages or 
PL/M programs . 

In general, the following steps must be observed 
to use the Floating-Point Library: 

1) An area of memory must be reserved for the Floating- 
Point Record (FPR) . 

(2) The names of the FPAL procedures must be declared to 
be "external". 

(3) FPAL procedure references must be imbedded in our 
source code where appropriate. 

(4) The FPAL procedure used by our program must be linked 
to our object file. 

(5) If we want to use the FPAL procedures in our program, 

that program cannot use symbols that are reserved for 
FPAL such as symbolic names beginning with a 
"commercial at" sign, or names whose second 

character is "Q" or "?". 

Floating Point Record (FPR) procedures, 
arithmetic procedures, error handling procedures, and inter- 
face to FPAL are described in Reference 2 in detail. 

Data Formats : 

FPAL procedures operate on single-precision 
binary numbers, either 32-bit integer format or in a 32-bit 
floating-point format. 
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The integer format recognized by the FPAL is a 
positive or negative (two's complement) 32— bit binary number 
The approximate range of this format is from -2.147x10^ to 

9 

+2.147x10 . Single-precision formats in the floating-point 



accximulator and 8080 memory areas are shown in Figure 6 




where s = sign bit 
e = exponent 

f = fraction or mantissa 
Figure 6. SINGLE PRECISION FORMAT FOR FPAL 



The three fields within these formats are: 

® ..is the sign bit. Sign magnitude repre- 

sentation wehre s = 0 means positive and 
s = 1 means negative. 

® is exponent bits. The exponent is offset by 

7 

(2 -1) . All zeros and all ones in the 
exponent field are currently reserved for the 
floating point zero and the invalid numbers 
described above. 
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f 



is fraction bits. When the exponent is nonzero 



a 1 bit is assumed at the left of the fraction; 



the binary point is between the assumed bit 
and the explicit fraction bit. 



The number base for the FPAL is binary. The 
value of a given binary representation can be formulated as 



operations, the correct fraction results but the exponent 
is "wrapped around". A "wrapped around" exponent is defined 
to be e^, where the true (offset) exponent, e^, can be 
derived from e^ by considering an expanded range of 
exponents and 




• (1. + .f) 



where e ^ 0 and e FF. 



When overflow or underflow occurs during FPAL 



on overflow 



e 



t 



e 



w 



( 3 . 2 ^ - 2 ) 



on underflow 



e 



t 



e 



w 



(3. 2^-2) . 



The floating point arithmetic library (FPAL) complete 
functions and execution times are given in Table 2. 
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FPAL Procedure 


Bytes 


Typical 

Execution 

Time 


Float. -Point Add. (FADD) 


463 


700 


Float. -Point Subtract (FSUB) 


463 


700 


Float. -Point Multiply (FMUL) 


404 


1500 


Float. -Point Divide (FDIV) 


342 


3600 


Conversion from decimal float. -point 
minber to binary (FQFD2B) 


■ 725 


not given 


conversion fron binary float. -point 
mitiber to decimal (FQFB2D) 


1585 


not given 


Float-to-Fixed Conversion (FIXSD) 


178 


not given 


Fixed- to- Float conversion (FLTDS) 


139 


not given 


Floating Point COrtpare (FCMPR) 


159 


300 


Floating Pint Test (FZTST) 


56 


not given 


floating Point Absolute (FABS) 


36 


not given 


floating Point Negative (FNBG) 


43 


not given 


Note: All time values are specified in 

microseconds . 



Table 2. 8080/8085 FPAL FUNCTIONS AND EXECUTION 

TIMES 
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c. 24-Bit Special Floating-Point Arithmetic 

Software Subroutines 

This software subroutine has been written in 
Reference 11 to perform the floating-point arithmetic. But 
the format that has been used is not our concern and not 
suitable to our standards. The only reason it was mentioned 
here is to show how slow execution times it has even with 
16-bit mantissa. Its functions and execution times are 
given in Table 3 . 

d. The SBC 310 High-Speed Mathematics Unit 

The SBC 310 High-Speed Mathematics Unit is a 
member of a complete line of Intel SBC 80 system expansion 
modules. In performing high-speed mathematic functions, 
the Math Unit acts as an intelligent processor slaved to 
one or more SBC 80 computer masters. The math unit performs 
its repertoire of 14 arithmetic functions an order of 
magnitude faster than is possible with software routines. 

Its functions and execution times are shown 
in Table 4 and a detailed description of programming infor- 
mation and principles of operation are explained in 
Appendix A because it was selected to perform this thesis. 

Four different systems and their features have 
been explained individually as they're chosen to be 
compared in the beginning of this section. The comparison 
among the SBC 310, High-Speed Math Unit, the AM9511 
Arithmetic Processing Unit, the Floating-Point Arithmetic 
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OPERATION NAME 


MAX EXECUTION TIME 


24-Bit Floating-Point Addition 
(ADD) 


1300 


24-Bit Floating-Point Subtraction 
(SUB) 


1400 


24-Bit Floating-Point 
Multiplication 

(MULT) 


3400 


24-Bit Floating-Point Division 
(DIV) 


4100 


24-Bit Floating-Point Comparison 
(COMPARE) 


800 


24-Bit Floating-Point Squareroot 
(SORT) 


17, 000 


24-Bit Floating-Point Cosine 
(TRIG 1) 


18,000 


24-Bit Floating-Point Sine 
(TRIG 2) 


18r 000 


24-Bit Floating-Point Cos and Sine 
(TRIG 3) 


30,000 


24-Bit Floating-Point Arctangent 
(TRIG 4) 


19,000 


Note: All time values are specified in microseconds. 



Table 3. 24-BIT SPECIAL FLOATING-POINT ARITHMETIC 
OPERATIONS NA^IE AND EXECUTION TIMES 
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OPERATION NAME 


OPERATION 

CODE 


TYPICAL 

EXECUTION 

TIME 


MAXIMCM 

EXECUTION 

TIME 


Fixed-Point Multiply (MUL) 


0 


15 


20 


Fixed-Point Divide (DIV) 


1 


26 


30 


Extended Fixed Point Div. (EDIV) 


E 


84 


100 


Float. -Point Multiply (FMUL) 


2 


84 


100 


Float. -Point Divide (FDIV) 


3 


92 


110 


Float. -Point Add (FADD) 


4 


33 


75 


Float. -Point Subtract (FSUB) 


5 


33 


75 


Float. -Point Square (FSQR) 


6 


84 


100 


Float. -Point Square Root (FSQRT) 


7 


178 


205 


Fixed-to-Float. Conversion (FLTEB) 


8 


72 


100 


Float. -to-Fixed Conversion (FIXSD) 


9 


42 


85 


Float. Point Catpare (FOIPR) 


A 


7 


7 


Float. Point Test (FZTST) 


B 


7 


7 


Exchange (EXCH) 


F 


4 


4 


Note; All time values are specified in microseconds. 




Listed times do not include time to pass arguments 


to the 


MATH UNIT and to read results tpon completion; this is 
typically 20 microseconds. 



Table 4. SBC-310 HIGH SPEED MATH UNIT FUNCTIONS AND 

EXECUTION TIMES 
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Library (FPAL) , and 24-bit Special Floating-Point Software 
execution times of their operations are shown in Table 5. 

The AM9511 APU has the fastest and SBC-310 High 
Speed Math Unit has the second fastest execution times. 
Also, if we compare the cost of the units, the AM9511 is 
only $195 in the market and it is cheaper than SBC-310 
Math Unit. But, when the purpose of this thesis concern, 
as a Kalman Filter, and immediate availability at the 
Naval Postgraduate School, the SBC 310 High-Speed Math Unit 
was chosen to perform the arithmetics of the Surface/ 
Subsurface Bearing-Only Tracking Radar algorithms. 

C. PHILOSOPHY OF SOFTWARE DESIGN 

The system was designed with the following objectives: 

1) Human engineered user interface. 

2) Capability of solving maneuvering board problems 
using Kalman Filter algorithms. 

3) Capability to display the result and the error 
performance of the assumed problem definition. 

In order to avhieve these objectives, the following 
modules were developed: 

a) "INT$DATA": This module contains range and bearing 

information of selected contacts and their observation time 
intervals for each scan. 

b. "MAIN$PR0G": This module is used to perform the 

implementation of Kalman filter algorithms, results and 
error calculations. 
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SYSTEM 

NAMES 

OPERATiai 

NAMES 


TYPICAL EXECUTION TIMES 


SBC 310 
High Speed 
Math Uhit 


AM9511 
Arithmetic 
Proc. Ifriit 


FDAL 

(Software) 


24-Bit 
Special F.P. 
Software 


32HDit Fixed Point Operations 


FMUL 


- 


97 


- 


- 


DIV 


84 


98 


- 


- 


32-bit Floating Point Operations 


EMUL 


84 


73 


1500 


3400 


FDIV 


92 


77 


3600 


4100 


FAIX) 


33 


27 


700 


1300 


FSUB 


33 


35 


700 


1400 


FSQR 


84 


- 


- 


- 


FSQRT 


178 


391 


- 


17,000 


Fixed-to-Float 

Conversion 


72 


28 


- 


- 


Float-to-Fixed 

Conversion 


42 


8 


- 


- 


FCEMPR 


7 


- 


300 


800 


SINE 


3800 


1898 


- 


18,000 


oosiEi: 


3800 


1920 


- 


18,000 


ARCTG 


5000 


2496 


- 


30,000 


Note: All time values are specified in microseconds. 



Table 5. COMPARISON OF FOUR DIFFERENT SYSTEMS EXECUTION 
TIME OF OPERATION 
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c. "DISPLAY"; This is used to display the result 
and error calculations on the CRT. 

d. "FLOATING$POINT" ; This module is used to perform 
all the necessary floating-point operations and to calculate 
the cosine and/or sine of a given angle in radians and the 
arctangent of the ratio of two input parameters. 

The whole program is listed in Appendix E. The formats, 
units and conversion factors used in the program are shown 
in Table 6. The Kalman filter will be described in 
Chapter IV. 
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FORMAT 


UNITS 


I/O 


INTERNAL 


I/O 


INTERNAL 


TIME 


XXXX . XX 


F.P. 


seconds 


same 


COURSE 


XXXX . XX 


F .P. 


degrees 


radians 


SPEED 


XXXX . XX 


F.P. 


knots 


n .mile/ 
sec . 


• 

BEARING 


XXXX . XX 


F.P. 


degrees 


radians 


RANGE 


XXXX . XX 


F.P. 


yards 


nautical 

miles 



X Numeric Character 

F.P Floating Point Represen-ation 



1 nautical mile 2025.3716 yards 

1 degree 0.0174532925 radians 

1 knot 1 nautical mile/hour 

7T 3.141593 



Table 6. FORMAT, UNIT AND COIJVERSION FACTORS USED 
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IV. KALMAN FILTER 



A two dimensional Kalman tracking filter provides a 
convenient framework for determining optimal filter 
parameters for x,y tracking for a two-dimensional radar 
measuring range R and bearing 0. 

A. DISCRETE KALMAN FILTER ALGORITHM 

The vehicle dynamics can be represented in matrix form 
as stated earlier. 



n+1 n n 

where Z , r and a were defined in Chapter II. 

The discrete form of a Kalman filter rather than 
continuous form is appropriate because of surface search 
radar antenna rotation results in scan while measurements 
of the contact position. The matrix form of the observation 
equation can be written as 



w 



n 



HZ„ + V 
n n 



(Ref. 13) 



where 



w 



n 



X (n) 
m 

y_ (n) 
^ m 
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X (n) = measured x coordinate at scan n 
m 

y (n) = measured y coordinate at scan n 
-‘m 



10 0 0 

H = 

0 0 10 



Z = Vehicle state vector at scan n 
n 



V 

n 



V (n) 

X 

V (n) 

y 



V (n) = random noise on x measurement at scan n, x (n) 

^ m 

V (n) = random noise on y measurement at scan n, y (n) . 

y m 



Optimal estimates of the contact state vector at 
scan n, are given by 



Z = Z + K (w - HZ ) (Reference 13) 

n n n n n 



where 

Z is the estimated state vector at scan n, 
n 

Z = $ Z , 

n n-1 

$ (see Chapter II, Section D) 
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= Optimal estimate of target state vector at scan (n-1) 



and 



n 



X 



n 



n 



n 



n 



n 



X 

r 

t 

X 



n 



n 



n 



n 



Z is the optimum estimate of the state vector after the 
n 

measurement w is processed 
n 



Z is the optimum estimate of the state vector before 
n 



the measurement is processed 



*• T ~ T -1 

K = P H (HP H + R ) 
n n r n 



(Reference 13) 



is the gain matrix at scan n and computed for steady- 
state as 



XX 



xy 



K 



n 



B /T 
XX 



B /T 
xy 



yx 

B /T 
yx 



yy 

B /T 

yy 
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P is the estimated covariance matrix of the estimation 
n 

errors prior to the processing run. 



R is the covariance matrix for the observation noise, 
n 



R 



n 



2 , , 



2 

a 

xy 



(n) 




2 

a 

y 



(n) 



where using the polar coordinate system shown in Figure 7. 



X = R sin 9 
m 



y = R cos 9 
■^m 



(WEST) 




(EAST) 
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The elements of the covariance matrix for observation 



noise are: 

a ^(n) = a ^ sin^ 0 + R^(n) cos^ 0 (n) 

X 3T_ y 

a ^(n) = a ^ cos^ 0 (n) + (n) sin^ 0 (n) a ^ 

y- V 0 

a (n) = ^ sin 20 (n ) [a ^ - R^ (n) 

xy 2 V / I j. \ / g j 

2 

= variance of the range measurement noise 

2 

a. = variance of the bearing measurement. 

0 

» 

R(n) and 0 (n) are the vehicle range and bearing at scan n. 
The covariance matrix can be computed recursively as 

m m 

= $ Pj^ <6 + rQr (Reference 13) 

with P = (I - K H) P 
n n n 

where P is the covariance matrix of estimation error 
n 

after processing w^. 

2 

Q = a and it is the variance of random contact 

3i 

acceleration. The variance of contact acceleration must 
be equal and independent along x and y axis, and from scan 
to scan. 
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This form of the Discrete Kalman filter is valid for 
our contact maneuver problem. But the solution of the 
equations in matrix form is unnecessarily complex for the 
microcomputer. Therefore, an algebraic solution of each 
element in the matrix form of equations was formulated. 
These equations are listed in Appendix D. 



B. KALMAN FILTER PARAMETERS 

The Kalman Filter performance is determined by three 

parameters : 

2 

1) ~ Variance of the range measurement error 
2 

2) Oq - Variance of the bearing measurement error 

2 

3) a - Variance of the contact acceleration. 

3 . 

The first two parameters are functions of the selected radar 
and the ship's gyrocom.pass . The third parameter is a 
function of the contact relative motion and the desired 
smoothness of the output data. 

In order to solve the tracking problem and figure out 
the parameters of the Discrete Kalman Filter, it is necessary 
to make measurements of 

a. Contact Range, 

b. Contact Bearing, 

c. Own Ship's Course, 

d. Own Ship's Speed, and 

e. Time of the observation. 
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It is assijmed that all noise is Gaussian and Received 
signal to noise ratios are large (>>1) . 

The minimvun performance of the selected radar (AN/SPS-10) 
characteristics are listed in Appendix C and the video 
signal to noise ratio with single pulse detection is given 
as 



S 

N 




(47t) 



"3 4 

kT F B L R 
o n n s 



where 



(Reference 8) 



= peak transmitted power (Watts) 

G = antenna gain 
X = wavelength (m) 

2 

a = radar cross section of contact (m ) 
a = attenuation constant of propagation medium (dB/m) 
R = range to contact (m) 

-23 

k = Boltzman's constant (1.23x10 joule/deg) 

= standard temperature (280°K) 

B^ = receiver noise bandwidth (Hz) 

Lg = system losses. 



The contact range is determined by measuring the time 
required for a radar pulse to travel to and from the contact. 



R = 



ct 

2 
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where c = speed of light 

t = elapsed time from pulse transmission to 
pulse detection 
R = range to contact. 

Threshold detection could be used to indicate the 
occurrence of the leading edge of the received radar video 
pulse. Leading edge detection is selected because it 
eliminates errors resulting from variations in the pulse 
width. The standard deviation of range error is given 
below if the rise time of the pulse is limited by receiver 
bandwidth B. 



R 



2B(2 |) 



1/2 



(Reference 8) 



The contact bearing is determined by measuring the 
bearing of the radar beam at the time the contact video 
is detected. Using beam splitting technique to estimate 
the center bearing of the contact, the standard deviation 
of bearing error is 

1.066 

= s ^T/2 

[N (|)^] 
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where 



a 



0 



Standard deviation of bearing error 




two way beamwidth 



signal to noise ratio at center of beam 



c 



N 



g 



number of pulses emitted as the antenna rotates 



20 



g* 



The standard deviation of range error and bearing error 
are functions of the video signal to noise ratio that is 
given on page 48 for a single pulse detection. When 
considering the threshold to noise ratio, it is necessary 
to define both the desired mean time between false alarms 
and minimum probability of detection of the contact. 
Probability of false alarms is given by Ref. 8 in terms of 
the mean time between false alarms and the IF bandwidth of 
the receiver as shown below. 



The mean time between false alarms is a function of 
the threshold to noise power ratio and the receiver IF 
bandwidth. And it can be increased if the tracking system 
is off whenever the contact is not in our selected range. 
The maximum unambiguous range of AN/SPS-10 Radar is 129.6 
nautical miles and for example 50 nm of that is of our 
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concern. This would disable the measurement circuitry 
61.4 percent of the time and it will give us an additional 
advanatage . 

The mean time between false alarms of 35.3 sec (10 times 

the mean interval) has been selected. The probability of 

—9 — 8 

false alarms 5.66x10 for one contact and 1.93x10 for 

50 nm of our interest range. And this leads to a threshold 

to noise power ratio of 12.5 dB. Also the reflected radar 

signals from contact will fluctuate due to cancellation 

and reinforcement of the waves reflected from various 

positions of the contact. These fluctuations can be formed 

as Swerling's Case I and additional 11 to 14 dB will be 

required to insure a probability of 0.9 to 0.99 resepectively . 

The standard deviation for bearing and range can be 
related to contact range and cross sectional area as given 
in Ref. 9. Based on a threshold to noise ratio of 12.5 dB 
and interested range of a 50 nm the standard deviation 
of the radar range and bearing errors will be 140.2 yards 
and 0.2062 (0.003599 rad), respectively. 

In addition to radar measurement errors also there are 
ship's gyrocompass and pitometer log errors which come from 
the own ship's course and speed measurements. And gyro- 
compass not only determines own ship course but also deter- 
mines the true bearing of the radar antenna. It is assumed 
that the standard deviation of own ship's gyrocompass 
error is 0.1 degrees (0.001745 rad) and standard deviation 
of own ship's pitometer log error is 0.1 knots. 



51 



Therefore, the variance of bearing error is equal to 

2 

the Siam of the radar bearing error (0.0425 deg ), and the 

2 

variance of the own ship's gyrocompass error (0.01 deg ) 

2 

leads to the total variance of bearing error 0.0525 deg 
“6 2 

or 15.998x10 rad . The variance of range error is 
2 

0.0048 nm . 

The variance of target acceleration parameter is the 
ability of the filter to provide the closest estimate of 
a contact's position and velocity for a moving contact. 

It is good to allow the experienced operator to choose the 
variance of acceleration parameter to be used with each 
contact. 

C. KALMAN FILTER INITIAL CONDITIONS 

As it is explained before the Kalman Filter algorithms 
represent the recursive calculations of the optimal 
estimates of the covariance of estimation error matrix and 
the contact state matrix. The initial condition of the 
contact state matrix is equal to the expected value of the 
initial state. The expected values of the initial position 
elements are the first measurements: 

Xq = R(0) sin 6 (0) 

Yq = R(0) cos 6(0) 

The expected values of the initial velocity elements are 
zero. 
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The initial condition of the covariance of estimation 



error matrix is equal to the covariance of error of the 
initial state estimate: 



P(0) = E[Zq - Z(0)] [Z(0) - Z(0)] 






P(0) = 



0 

0 

0 



E[(Xq) 

0 

0 



0 

0 



E[(Yq-Yq) ] 



0 

0 

0 



E[Yq)^] 



The time required for the Kalman Filter to achieve 
a steady-state condition is very sensitive choosing the 
closest right values of the nonzero elements of the initial 
state estimate matrix. They were chosen to be brought the 
contact on y-axis as an initially by given equations in 
Reference 13. 



53 



V. SYSTEM DESCRIPTION 



The description of the system is divided into three 
major areas: hardware dependencies, system characteristics, 

and system performance. 

A. HARDWARE DEPENDENCIES 

Because of hardware equipment that was selected the 
following hardware dependencies exist in the current 
implementation of the system: 

1. The system utilizes an SBC-310 High-Speed Mathematics 
Unit to perform floating point arithmetic. Although as 
discussed before the presence of this math unit could be 
avoided by replacing its functions with appropriate software 
routines performing the same operations using the same 
formats, this is not recommended because of excessive over- 
head that would result, especially with regard to the 
execution time. Appendix F explains how the Math Unit was 
actually implemented. 

2 . The system depends in three ways on the type of 
terminal used. The serial asynchronous procedure necessary 
to communicate between CPU and the terminal, the code needed 
to control the CRT's functions, and the general features 

of the DATAMEDIA Elite 2500 Video Terminal notably the 
programmable roll mode, the setting of privileged fields, 
the capability of having an addressable cursor, and the 
possibility of making displayed messages blink. 
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3. The system occupies approximately 13K bytes of 
physical memory for code, and approximately 3K bytes to 
be used for variable data, therefore a configuration of 
at least 16K bytes of RAM is necessary to execute the 
system. 

B. SYSTEM CHARACTERISTICS 

As it was established before, the system was designed 
to perform basically to solve the maneuvering board problem 
finding the contact speed and course. 

Due to interaction capability between modules as 
allowed by the language PL/M-80 through the use of attri- 
butes PUBLIC and EXTERNAL for the procedures, the following 
list was written in order to show interactions. This list 
shows the modules which have procedures called by the ' 
listed module. 



1 . EXTER 

2 . INT$DATA 

3. MAIN$PROG 

4 . DISPLAY 



a. 


NUMOUT 


b. 


DISPLAYS CRT 


c. 


CRTS READ 


d. 


checkSyesSno 


e. 


DISPSERROR 


FLOATINGS POINT 


a. 


MUL 


b. 


DIV 


c . 


EDIV 


d. 


FMUL 


e. 


FDIV 


f . 


FADD 


g- 


FSUB 


h. 


FSQR 


i . 


FLTDS 


j • 


FIXSD 
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k. FSQRT 

l . FCMPR 

m. FZTST 

n. COS$SIN 

o . ARC$TAN 

The Module "INT$DATA" is used to implement the radar 
measurements: range as a yard, bearing as a degree, time 

of measurements as a second and their interfaces process 
the system. 

The "MAIN$PROG” Module used to perform the Kalman filter 
algorithms, computation of the filtered target course and 
speed, display of the result and error performance. 

The Module "DISPLAY" includes' ithe display procedures 
on the CRT 

The "FLOATING$POINT" Module is used to perform the 
floating point arithmetic. 

In order to achieve these objectives, a bottom-up 
implementation philosophy was chosen. Because of the 
modular design encouraged by PL/M 80 and ISIS-II, it was 
also possible to map the different levels of design into 
corresponding modules of software. The basic idea was to 
encompass all functions corresponding to a level of design 
into one software module capable of performing all the 
necessary functions. 

C . SYSTEM PERFORMANCE 

The performance of the Kalman filter was evaluated for 
maneuvering contact. The own ship was proceeding on a 
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steady course (000°) North and at a fixed speed of 30 knots 
from an initial position at the origin of the cartesian 
coordinate system as shown in Chapter II, Figure 3. The 
contact was heading to 060° true (North/East) by 36 knots 
and its initial position was 009.6° (North/East), 37,800 
years from own ship. 

The scenario has been made to change the contact 
course and speed after every 15 sets of measurements. The 
approximate problem was solved on the maneuvering board and 
by hand calculator. The first chance has been made for 
contact course and speed which were 073° (North/South) true 
and 50 knots, respectively, then second reorientation has 
been made to 328° (North/West) true course and 30 knots 
speed. 

The accuracy of the computed course and speed of the 

contact will depend on the tracking system parameters and 

the maneuverability of the contact. The antenna rotation 

rate determines the mean time between measurements, which 

in the case of the AS-1161/SPS Antenna equates to 3.53 

seconds. The majority of surface contacts will be capable 

of turning radii of 100-1000 yards, velocities of 0-50 

2 

knots, and acceleration of 0-0.5 yards/sec . The following 

2 2 

Kalman Filter measurement parameters (a. and a ) were 
used as calculated in Chapter IV, Section B. 



2 



0.0048 nm 



2 



a 



r 



2 



15.988x10 



-6 



rad 



2 



a 



e 
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The Kalman filter initial conditions were the same as 



those discussed in Chapter IV, Section C. The variance 

2 

of contact acceleration (a ) parameter was selected as 

2 2 4 

a = 0.0004 nm /sec 

5L 

Also as is mentioned before, the variance of contact 

2 

acceleration (a_ ) parameter is an independent variable and 

a 

it is a compromise between the user's desired smoothness 
of the output data and the ability of the filter to provide 
the closest estimate of a contact's position and velocity. 

This parameter could be selected by enabling one of 
several read-only memories (ROM) containing the desired 
values. The particular ROM installed at any one time. could 
be selected from a library of such ROM's by the Commanding 
Officer, depending on the ship's mission. 

The filtered outputs of the contact course and speed 
according to the selected scenario were plotted on Figure 8 
and Fiture 9, respectively. As is noted on those figures, 
the Kalman Filter achieved a steady state condition after 
four to seven contact position measurements, corresponding 
to 14.1 and 24.7 seconds related to the amount of reorienta- 
tion has been made on contact course and speed. Steady state 
is defined as the condition where the Kalman Filter Gain 

Matrix (K ) would be constant for contact with constant 
n 

relative position and constant sample intervals. 
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Figure 8. CONTACT COURSE VS. TIME 
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Figure 9. CONTACT SPEED VS. TIME 



The bearing, range, course, and speed of contact was 
computed from the Kalman filter's optimal estimate of 
contact state and compared with the actual contact bearing, 
range, course and speed. Figures 10, 11, 12, and 13 show 
plots of the filtered output of the contact bearing, range, 
course and speed errors, respectively, as a function of 
the time interval between measurements. 

Several conclusions can be made from these figures. 

First of all the Kalman filter initial conditions are 
very much affected on bringing the Kalman filter to steady 
state condition. Any information about contact speed would 
cause to shorten the time spent achieving a steady state 
condition. Also selecting the variance of target accelera- 
tion parameter is as important as initiating the filter. 

The value of the variance acceleration parameter 0.0004 
2 2 

nm /sec were chosen after many trials of that selected 
scenario. It should be noted that the error averages for 
contact are a function of the elapsed time between maneuvers 
and the time spent in maneuver . 
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Figure 10. BEARING ERROR VS. TIME 
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Figure 11. RANGE ERROR VS. TIME 
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Figure 12. COURSE ERROR VS. TIME 



7.5 




(sq.oui{) aojoa peads • 



65 



Figure 13. SPEED ERROR VS. TIME 



VI . CONCLUSION 



A microcomputer-based tracking system has been developed 
that eliminates an inaccurate and tedious task of manual 
tracking, for any ship, without the expense of NTOS size 
aquipments, with digital automatic surface/subsurface capa- 
bility. 

It is possible to obtain optimal estimates of contact 
position and velocity, with the aid of the 4-state Discrete 
Kalman Filter and appropriate conversions for polar coordinate 
neasurements . Thus the contact course and speed computations 
vere successfully performed by using the Intel intellec MOS 
nicrocomputer system, much more accurately than manual 
tracking techniques. 

The SBC-310 High Speed Mathematics Unit has been selected 
based on a standard floating-point number format with 24 
Dit mantissa, 7 bit exponent and a sign bit by comparing 
sxecution times of procedures to perform sufficient accuracy 
to reliably compute the Kalman Filter Algorithms and con- 
(/ersion factors. 

The total average of the system execution time was 
approximately 50 msec and with the 3.5 sec sweep of an AS~ 
1161/SPS radar antenna, a single system could continuously 
provide bearing, range, course and speed data for more than 
fifty contacts by modifying the software program. 
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For future work it is possible that the system could 
compute a contact's closest point of approach (CPA) and could 
display as (1) CPA Bearing, (2) CPA Range, and (3) Time of 
CPA by a more sophisticated software study. 

The system can be a valuable tool aboard non-NTOS ships 
as implemented. A microprocessor-based system such as the 
one described, has the potential to provide any ship with a 
digital surface/subsurface tracking capability without the 
expense of NTOS size equipment. It can reduce the manning 
of underway watches while improving the quality of the 
tactical information available to the Officer of the Deck. 
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APPENDIX A 



Intel SBC-310 High Speed 
Mathematical Unit Descriptions 

A. GENERAL INFORMATION 

The floating-point package developed for this system is 
based on the SBC-310 High-Speed Mathematics Unit fromintel 
Corporation. As described by Reference 1, the SBC 310 Unit 
is a member of a complete line of the Intel SBC 80 System 
expansion modules. In performing high-speed mathematical 
functions, the Math Unit acts as an intelligent processor 
slaved to one or more SBC 80 Computer Masters. The Mathe- 
matic Unit performs its repertoire of 14 arithmetic functions 
an order of magnitude faster than is possible with software 
routines . 

B. DESCRIPTION OF THE MATH UNIT 

The Math Unit is a microprogrammed processor on a single 
board and is designed to be plugged into a standard SBC 
604/614 Modular Backplane and cascade to interface directly 
with an SBC 80 single board computer or to be used with an 
Intel Intellec Microcomputer Development System CMDS) . 

The Math Unit includes the following standart Intel Series 
300 shotthy bipolar components: 3001 Microprogram Control 

Unit (MCU ) , 3002 Central Processing Element CCPE) , 3003 Look- 
Ahead Carry generator CLCG) , and 3604 Electrically Programmable 
Read Only Memory CPROM) . Also included are pipe line register 
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and bus interface logic. The pipeline register permits the 
overlapping of microinstruction fetch/execute cycles and the 
bus interface logic provides compatibility with the Intel 
Multibus . 

Standard Operations include floating point add, subtract, 
multiply, divide square and square root; fixed point integer 
multiply, divide, and extended divide; conversion between 
fixed and floating point representations; and test, compare, 
and argviment exchange operations . 

The Math Unit implements unbiased rounding for maximum 
accuracy. Unbiased rounding is the same as ordinary rounding 
unless the result is exactly midway between two floating 
point nxambers; in this case ordinary rounding always increases 
the result, whereas unbiased rounding rounds the result to 
the nearest even number. When a calculation is performed 
that results in either an exponent underflow or overflow, the 
Math Unit provides exponent wraparound to prevent loss of 
information. 

Operation Codes for invoking the arithmetic functions are 
passed to the Math Unit via I/O Write, Commands, which are 
also used to initialize the unit with a memory base address. 

I/O read commands are used to determine the math unit status. 
Arguments are passed to the Math Unit via Memory Write, 

Commands and the results are obtained via Memory Read Commands, 

The Math Unit which can be operated either in the Inter- 
rupt of Polled mode, generates a busy signal during processing 
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operations and generates either a complete signal or an Error 
signal after the computation is complete. The information 
to the host computer which these three signals convey is 
explained in Ref. (1) . 

The memory base address and I/O base address are user 
selectable. The 16-bit memory address is completely under 
software control and is assigned by the host processor through 
a sequence of I/O Write Commands, addressed to the Math Unit. 
The 8 bit I/O base address is selected by a dual inline 
package (DIP) switch on the board. 

All Math Unit operations including arithmetic calculations, 
data flow between functional elements on the board bus inter- 
face, and associated logical tasks, are resident microprogram 
permanently stored in a set of eight Intel 3604 Erasable 
Programmable Read Only Memory CEPROM) chips. This memory 
provides 1,024 micro-instructions of 32 bits each. 

C. PREPARTION FOR USE 

1 . Pnstallation Consideration 

The Math Unit is designed for interface with an Intel 
SBC 80 single board computer based system or an Intel Intellec 
Microcomputer Development System (MDS) . 

When installing the SBC 310 in an Intel Intellec MDS, 
the CPU board needs to be reconfigured in order to generate 
a Qualified Write Signal; this reconfiguration was obtained 
by considering the advanced acknowledge (AACK) feature by 
moving a jumper labelled "advanced write" from a D-C connection 
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to an E-D connection and by disabling the AACK/LINE (pin 25) 
on the CPU board as shown in the schematic diagram on page 
3-47 of Reference 4. 

2 . I/O Base Address Switches 

The host processor transmits control information 
and receives status information from the Math Unit by issuing 
I/O Write and I/O Read Commands, respectively. The I/O 
address used for these commands is relative to an 8-bit base 
address that must be a multiple of 8. This base address is 
assigned by the user by means of an 8-pole dual inline package 
(DIP) switch assembly. Five of the eight switch poles are 
ronnected to the I/O base address detection logic; the other 
:hree poles are unused. 

The Math Unit had used its DIP switch set to the 
C/0 base address of 10 hexadecimal (only switch pole no. 4 
<;as set on) . 

3 . Programming Information 

The I/O base address, which must be assigned by switch 
selection before the memory base address can be assigned, is 
lormally performed as part of the initial installation proce- 
iure. This switch setting allows the user to establish a 
reference or base to the ports being used in the I/O opera- 
tions; Table 7 shows the configuration of the I/O addressing 
as it was set in the system. 

The memory base address, which is software controlled, 
is assigned by a sequence of two I/O Write Commands. The 



71 



I/O PORT 
ADDRESS 


OUTPUT 


INPUT 


P = OlOH 


OP CODE 


R 


P+1 = OllH 


MEM. LOW 


STATUS BYTE 


P+2 = 012H 


MEM. HIGH 


R 


P+3 = 013H 


R 


R 


P+4 = 014H 


R 


R 


P+5 = 015H 


R 


R 


P+6 = 016H 


R 


R 


P+7 = 01 7H 


R 


FLAG BYTE 


P: I/O BASE ADDRESS 
R: RESERVED 

OP CODE: MATHEMATIC FUNCTION; See Table 4 
MEM. LOW: Memory BASE ADDRESS (Lower Byte) 
MEM. HIGH: Memory BASE ADDRESS CUpper Byte) 



Table 7. I/O ADDRESSING OF SBC-310 HIGH SPEED MATH UNIT 
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first command is addressed to port (P+1) and loads the low 
order byte of the memory base address. The second command is 
addressed to port (P+2) and loads the high order byte of the 
memory base address. The memory base address must be a 
multiple of 16; i.e., the lower byte must be in the form XOH 
(X is any hexadecimal digit) to accommodate the 16 required 
memory locations used in the arithmetic operations . In the 
system developed for this thesis the memory base address was 
set of 0F790H. After both bytes are output, the memory base 
address (M) is established and need not be reloaded during 
any subsequent operations. An initialization routine for 
establishing the memory base address was designed and it can 
be seen in the "INIT$FP" procedure in the floating point 
module (see "FLOATING$PINT'' module in Appendix F) . 

4 . Math Unit Functions 

The Math Unit performs floating point arithmetic, 
fixed-point integer arithmetic, compare and test operations, 
and, float-to-fix and fix-to-float conversions. Operation 
Codes and execution times for the various functions are 
listed in Table 4. Arithmetic and conversion formats are 
shown in Table 8. 

Beyond the functions performed by the Math Unit, 
the floating-point package was designed with two more proce- 
dures which compute the cosine and sine of a given angle and 
the arctangent value of the ratio of given arguments . 



73 



SINGLE PRECISION FLOATING POINT 



Mf3 m-2 
S7 0 22 16 



M 

15 8 7 0 



=^9" -exponent fraction 

8 bits 23 bits 

where: M: Memory Base Address 

S: ”0" = positive; "1" = negative 

E7-E0: Biased exponent; Bias = 07FH. 

F22-PO: Fracrtion; F is always normalized 



FIXED POINT INTEGER 



M+1 

5 U- 



M 



binary point 
vdiere: M: Memory Base Address 

F15-F0: 16-bit integer (unsigned) 



M+3 Mf2 

31 2423 16 



Mfl 



M 






8 7 



>A 



where: M: Menory Base Address 

F31-F0: 32-bit integer (unsigned) 



binary point 



CONVERSION FUNCTION 



Mf3 M+2 
31 2423 16 



M+1 

15 8 7 



M 



Sign 

viiere: M: Memory Base Address 

S: "0" = positive; "1" = negative 

F30-F0: to^'s conplonent integer 



binary point 



Table 8. ARITHMETIC AND CONVERSION FORMATS FOR 
SBC- 10 HIGH SPEED MATH UNIT 
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5 . Argument and Result Data Formats 



Argument and result data formats and memory locations 
for the various operations are presented in Table 9. For each 
argument and result, this table includes a Format number 
cross-referenced to one of the four formats shown in Table 8, 
Table 9 also includes the OP CODE for each operation. It is 
important to note that the result of an operation replaces 
the first argument in memory, and that the second argument 
may be destroyed in the course of the computation, these 
side effects were avoided in the floating-point software 
design by saving the original values and by allowing the 
user the possibility of having one of the operands as the 
result. Error conditions for each operation are described 
in the next paragraph. 

6. Status and Flags 

The Math Unit may be operated in the interrupt mode 
or piled mode. 

In the interrupt mode, the Math Unit may be wire- 
wrapped to initiate an interrupt request under one or both 
of the following conditions: 

a. Operation Complete without an error. 

b. Operation Complete with an error. 

These "completion" signals may be individually wire-!- 
wrapped to separate interrupt lines or both "completion" 
signals may be wire-wrapped to the same interrupt line (Reference 
1 ) . 
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OP 


OP 

CODE 


ARGUMEOT 

FORMAT* 


ARGCMENTS** 


RESULT 

FORMAT* 


RESULT*** 


MUL 


0 


2 

2 


M,Mfl 

Mf4,Mf5 


2A 


M,Mfl,Mf2,Mf3 


DIV 


1 


1 

1 


M,Mfl 

Mf4,M+5 


2 


M,Mfl (ronainder 
in Mf4,Mf5) 


EDIV 


E 


2A 

2 


M,Mfl,Mf2,Mf3 

Mf4,Mf5 


2 


M,Mfl,Mf2,Mf3 
(ranainder in Mfl, 
M+5,Mf6,M+7) 


FMUL 

FDIV 


2 

3 


1 


M,Mfl,Mf2,Mf3 

Mf4,I4+5,Mf6,Mf7 


1 


M,Mfl,Mf2,Mf3 


FADD 


4 










FSUB 


5 










FSQR 


6 


1 


M,Mfl,^B-2,^^f3 


1 


M,Mfl,Mf2,Mf3 


FSQKT 


7 










FLTDS 


8 


3 


M,M+l,Mf2,M+3 


1 


M,Mtl,M+2,Mf3 


FIXSD 


9 


1 


M,M+l,Mf2,M+3 


3 


M,Mfl,Mf2,M+3 


FCMPR 


A 


1 


M,Mfl,M+2,Mf3 

Mf4,Mf5,Mf6,M+7 


- 


STATUS Byte 


FZTST 


B 


1 


M,Mfl,Mf2,Mf3 


- 


STATUS Byte 


EXCH 


F 


Any 


M,MH,Mf2,Mf3 

M+4,Mf5,Mf6,jyH-7 




Rotates both 
Arguments 



*Refer to appropriate FORMAT No oolunn in Table 8 



**Second argument is always the operator, and may be destroyed 
during operation. 

***Results of all operations, except FIXSD are rotinded, FIXSD 
truncates the result. 



Table 9. OPERATION ARGUMENT AND RESULT DATA FORMATS FOR 
SBC-10 
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I 
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I 




In the polled mode, the subroutine designed fo ithis 
purpose checks both the status byte and the flag byte (Tables 
10 and 11) . The polled mode procedure loops on testing the 
busy bit until the busy bit is clear and then checks the 
error bits. If an error exists, the error code is input from 
the Math Unit and a message error is issued. The software 
has developed for control of the Math Unit if the system 
employs the polled mode described. 

As mentioned before, the condition of the Math Unit 
is continuously updated and stored. The flag byte shown in 
Table 10 may be obtained by performing an I/O Read Command 
to Pt 7 (Table 7) . After an operation is completed, the 
status byte may be obtained by performing an I/O Read Command 
to P+L (Table 7) . 

As shown in Table (11) , the status byte indicates 
error conditions where applicable and the results of compare 
(FCMPR procedure) and Test (FZTST procedure) operations. 

Each of the six error conditions are defined as follows; 

a. Divide by Zero: (001) 

This error condition is returned by either "DIV" , 
"EDIV" , or "FDIV” procedures to indicate that an attempt was 
made to divide by zero. 

b. Domain Error (010) 

This error condition is returned by the "FSQRT" 
procedure to indicate that the argument was not in the domain 
of the function; i.e., an attempt was made to take the square 
root of a negative number. 
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7 


6 


5 


4 


3 


2 


1 


0 


R 


R 


R 


R 


R 


E 


C 


B* 



where : 

R is reserved for future use 
B is busy 

C is operation complete without error 
E is operation complete with error 



when B = 1, the Math Unit is busy and can not 
respond to further request except request for 
flags . 



Table 10, FLAG BYTE FORMAT FOR SBC-310 



7654 3210 





> 


< 


R 


R 


ERR 



where ; 

R is reserved for future use 
= is equal (for FCMPR and FZTST) 

> is greater than (for FCMPR and FZTST), 
< is less than (.for FCMPR and FZTST) 



ERR is a 3-bit error code specifying one of 
the following error conditions » 

QQO No error 
QQl Divide by Zero 
QIO Domain Error 
Oil Over flow 

100 Under flow 

101 First argument invalid 

110 Second argument invalid 

111 Reserved 



Table 11. STATUS BYTE FORMAT FOR SBC-310 



79 



c. Overflow (Oil) 

This error condition is returned by the "FADD", 
"FSUB", "FMIL", "FDIV", "FSOR", and "FIXSD" procedures. In 
the case of "FIXSD" procedure, this error indicates that the 
floating-point number is too large to be converted to a 32 
bit two's complement signed integer. If an overflow error 
occurs during "FIXSD", the floating-point argiament is left 
unchaged and may be read from the Math Unit. 

In all other cases, this error condition signifies 
that the exponent of the result is too large to be repre- 
sented in eight bits. In this case, OBEH is subtracted from 
the resulting exponent (bringing it back into range for other 
computations and ensuring a valid result) , and the lower 
eight bits of the exponent are returned in the exponent field 
of the result. 

d. Underflow (100) 

This error condition is returned by "FAOO", 

"FSUB", "FMUL", "FDIV", and "FSQR" procedures to indicate 
that the exponent of the result is too small to be represented 
in eight bits. In this case OBEH is added to the resulting 
exponent bring it back into range for other computations and 
ensuring a valid result, and the lower eight bits of the 
exponent are returned in the exponent field of the result. 

e. First Argument Invalid (101) 

This error condition is returned by the "FADD" , 
"FSUB", "FMUL", "FDIV", "FSQR", "FSQRT", "FIXSD", "FCMPR", 
and "FZTST" procedures to indicate that the first (or only) 
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argument for the specified function is invalid. The second 
argument (if applicable) is not checked if this error is 
encountered. The invalid argument is left unchanged and 
may be read from the Math Unit. 

f. Second Argxament Invalid (110) 

This error condition is returned by the "FDD", 
"FSUB", "FMUL", "FDIV" and "FCMPR" procedures to indicate 
that the second argument for the specified function is 
invalid. This error condition occurs only after the first 
argument is checked and found valid. The invalid argument 
is left unchanged and may be read from the Math Unit. 

Notes : 

(1) The floating-point argument may be expressed 

as : 

(1) X (2^ - BIAS) X (l.F) 

Notice that a "1" is assumed in the highest position of 
fraction. 

(2) There is one unique representation for zero; 

S = 0 

E7 - EO = 0 

F22 - FO = 0 
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(3) The following representations are invalid. 

(a) Assures unique representation of zero 

E7 - EO = 0 , and 

S 0 or F22 - FO = 0 

(b) Reserved for future enhancements 
E7 - EO = OFFH 

7 . Examples of Floating-Point Number~" Representations 

As shown in Table 8 the representation of a floating- 
point number has one particularity that needs to be men- 
tioned; i.e. , a ”1” is always assumed in the highest bit 
position, and then it yields an effective 24-bit mantissa. 

For the sake of clarity some examples are given in Table 12. 
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F.P. Number 


M 


M+1 


M+2 


M+3 


- INFINITY 


.OFFH 


OFFH 


OFFH 


OFFH 


- 255.0 


OOH 


OOH 


07FH 


OC3H 


5.0 


OOH 


OOH 


OAOH 


OCOH 


1 

• 

O 


OOH 


OOH 


080H 


08FH 


O 

« 

o 


OOH 


OOH 


OOH 


OOH 


+ 1.0 


OOH 


OOH 


080H 


03FH 


+ 5.0 


OOH 


OOH 


OAOH 


O40H 


+ 255.0 


OOH 


OOH 


07FH 


043H 


+ INFINITY 


OFFH 


OFFH 


OFFH 


07FH 


7T 


ODBH 


OOFH 


049H 


040H 


it 


OFFH 


OFFH 


07FH 


07FH 


itit 


OOH 


OOH 


08 OH 


OOH 


1.07 X 10^ *** 


OFFH 


OFFH 


07FH 


04BH 



* 

This is the largest number in the single-precision 
floating-point format, 

** This is the smallest positive number in the single 
precision floating-point format. 

*** This is the largest number that can be converted 
to floating point without losing accuracy. 



Table 12, FLOATING POINT NUMBERS 
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APPENDIX B 



System Start'-Up Procedure 

Caution: Never turn on or off the diskette drive with 

a diskette inserted I ! ! 

1. TURN ON MDS SYSTEM. Use the key located at the upper 
left corner of the front panel and turn it clockwise. 

The power indicator should light. 

2. TURN ON DISKETTE DRIVE. Use power switch located at the 
front panel. The on indicator should light. 

3. TURN ON DATAMEDIA TERMINAL (.CRT). Use switch located 
on right side. The cursor should appear at the screen 
after a few seconds. Ensure that the lights CD, CTS, 

ROLL and FULL DUPALLEX are on. 

4. PLACE SYSTEM DISKETTE IN DIRVE 0 and the source diskette 
in drive 1, with the read/write access slot first. Close 
door of the drives after insertion. 

5. BOOTSTRAP THE ISIS-II OPERATING SYSTEM. 

a. Press top of the Intellec BODY Switch. 

b. Press top of the RESET Switch. 

c. Observe that INTERRUPT 2 INDICATOR goes on before 
preceding. 

d. Press space-bar of Datamedia video terminal keyboard. 

e. Observe that INTERRUPT 2 indicator goes off before 
preceding. 

f. Press bottom of BOOT Switch. 
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g. Observe that the following message appears at the 
Datamedia Video terminal screen: 

ISIS-II, V3.4 

6. Issue the following Command; 

- ;F1; FATIH 

7. After a few seconds, the Datamedia Video terminal screen 
will show the first computation of COURSE, BEARING, 

RANGE, SPEED, THE KALMAN FILTER GAIN MATRIX ELEMENTS, AXX, 
AXY, AYX, AYY, BXX, BXT, BTX, BYY and the FILTER Error 
Performance DELTA RANGE, DELTA BEARING, DELTA RANGE DOT, 
DELTA BEARING DOT and following message appears at the 
screen ; 

NEXT ? CY/N) 

3. PRESS THE LETTER 'Y' and see ^YES*^ message appear on the 
screen to compute next measurement. 



85 



APPENDIX C 



AN/SPS-10 Search Radar Characteristics 
with AS-1161/SPS Radar Antenna 



Peak Power, 


190 KW 


Antenna Gain, 


30 dB 


Pulse Repetition Rate 


625 Hz 


Frequency 


5825 MHz 


Propagation Attenuation 


0.6 X 10 ^ dB/m 


Noise Figure 


22.9 dB 


Band width 


5 MHz 


Noise Bandwidth 


(1.05 X B) t 5.25 MHz 


System tosses 


5 dB 


Antenna Rotation Rate 


17 rpm 



Antenna Horizontal Beamwidth 1.9® 
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APPENDIX D 



Algebraic Form of Discrete Kalman Filter 

1. COVARIANCE MATRIX FOR THE OBSERVATION NOISE 

R 

n 

(n) = • sin^6(n) + r^ (n) • cos^SCn) 

a^(n) = • cos^6(n) + r^ (n) * a? • sin^6(n) 

y r o 

a^y(n) = j • sin 26 (n) • [a^ - r^(n)ag] 

2. THE GAIN MATRIX 

K = P(n) * H^ •(H-P(n)*H'^ + R )“^ 
n n 

A = (Pj^j^Cn) +a^(n) ) (P 33 (n) +ay(n)) 

- (P3i(n) + a^y(n)) (P3^3(n) + a^y(n)) 

A^x(n) = [Pj^j^Cn) (P33(n)+ay(n))-P^3(n) CP33^(n)+axy(n))J/A 
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II 

>1 

X 

< 


[Pj^3(n) +a^(n)) -Pn(n) (P13 Cn) +a^y(n))]/A 


A (n) = 

yx 


[P^j^ (n) (P33 (n) + Gy (n) ) - P33 (n) (P^^^ (n) + (n) ) ] /A 


> 

II 


[P33 (n) {P^j^ (n) + a^(n) ) - P^^ (n) (P33 (n) + (n) ) ]/A 




= 1^21 ^^33 (n)+Oy (n) ) -P23 (n) (P^j^ (n) +cr^y (n) ) ]/A 


B^(n)/T 


= tP23(n) (P3^j_(n)+a^(n) )-P23^(n) (P^_3 (n) (n) ) ] /A 


By^(n)/T 


= [P^^ (n) (P33 (n) +Oy (n) ) -P^3 (n) (P3^ (n) (n) ] /A 


B (n)/T 

yy 


= [P^3(n) (P^j^(n)+o^(n) )-P^j^(n) (Pj^3 (n) (n) ) ] /A 



3. THE COVARIANCE MATRIX OF ESTIMATION ERROR AFTER 
PROCESSING THE OBSERVATION 





P^ = U-K„H) P^ 




^11 <"><1 -\x> - ^31 f"' V 


^12 ” 


P 32 (n) Cl - -Pjjtn) Ajjy 


Pl3(n) = 


- ^33<"> V 


Pl^Cn) = 


^14‘"X^-*xx' - ^34^''’ V 



^ *XJ> 



Pjj^(n) " '’ll'"’ ®xx ■ '’ 31 *"’ ®xy 

22 (n) = P22(n) - Pj2 (n) - P32 (n) 
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4. THE COVARIANCE MATRIX OF THE ESTIMATION ERRORS 
PRIOR TO PROCESSING THE OBSERVATION 
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Ini tif p: 

procedure external* 
end ; 
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cos$sin : 

procedure (a,b,c) external* 
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ISIS-II PL/M-80 V3.1 COMPILATION OF MODULE INTDATA 
OBJECT MODULE PLACED IN :Fl :FATIH3.0BJ 
COMPILER INVOKED BY: PLM80 :F1 :FATIH3 .PLM 
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3. BXX/T = (PB21 * PB33Y - PB23 * PP31XY) / DELTA 

130 2 CALL FMUL ( . PB21 , .PB33Y , . BXX ) ; 

131 2 C<\LL EMUL ( . PB23 , .PB31XY , . TEMP ) I 
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/* 8. BYY/T = (PB43 * PBllX - PB41 ’S' PB13XY) / DELTA 

150 2 CALL EMUL ( . PB43 , . PBl IX , .BYY ) ; 

151 2 CALL EMUL ( . PB41 , . PB13XY , .TEMP ) ; 

152 2 CALL ESUB (. BYY , .TEMP , .BYY ) ; 

153 2 CALL EDIV ( .BYY , .DELTA , .BYY ) ; 
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164 2 C*LL FSUB ( . ? A13, . TEMP , . P A13) J 
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180 2 CALL EMUL f .PB14, .BXX , .PA24 ) ; 

181 2 CALL FMUL ^ . PB34 , . BXY , .TEMP ) ; 

182 2 CALL FADD ( . PA24 , . TEMP , . PA24 ) ; 

183 2 CALL FSUB ( .PB24 . .PA24 . . PA24) ; 



/* 9. PA:^1 = PB31 ONE^AYY - PBll AYX 

184 2 CALL EMUL ( .PB31, .ONE$AYY , . PA31 ) ; 

185 2 CALL EMUL { .PBl 1 AYX , .TEMP) ; 

186 2 CALL FSUB ( . PA31 , . TEMP , . PA31 ) ; 
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PL/M-60 COMPILER 
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/* 4. PB14 = P.»14 + PAP4 T 

228 2 CALL FMUL ( .PA24 , .T, .PB14) J 

229 2 CALL PAPD ( .PB14 , . PA 14 , . PB14 ) 
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/* 10. ?B32 = PA32 + Pft42 T 

243 2 CALL PMUL ( .PA42 , .T, .PB32) ; 

244 2 C«LL P*PD ( .PB32 , . PA 32 , . PB32 ) ; 
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/* 16. PB44 = PA44 + SIGAT2 
263 2 CALL PAT)D ( . PA 44 , . S IG AT2 , . PB44 ) ? 



I 





OPTIMUM PSTIMATK OF THF STATF VECTOR BEFORE THE MEASUREMENT. 
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/❖ VYFST* = VYESTB + BYY * DFLTAY BYX * DELTAX 

280 2 CALL FMUL ( .BYY , .DELTAY, . VYESTA) i 

281 2 C«LL FMUL ^ .BYX , .DELTAX, .TEMP) ; 

282 2 CALL FADD ( .VYESTA ,. TEMP , .VYESTA ) ; 

283 2 CALL FADD ( .VYESTA ,. VYESTB. .VYESTA) ; 
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318 2 CALL S^NDiCT^LF* 

319 2 CALL CHT$PRINT$STRING(.MSG$4); 

32t) 2 CALL DISPLAY$CRT( .SPEED); 

321 2 CALL SEND^CPLE; 







a. 




\ 



C!> 

2 

r: 















\ 




4^ 
















» 




















C3 




















P 




d5 
















P 


• •> 


p 
















z: 




p 
















CO 


o 


2 
















O 


p 


CO 
















o 


p 


O 


















w 


o 


















«r: 




• •»> 
















Eh 


45* 
















X c!3 O O 


P 




O O 














<c crj P p 


P 


>- 


2 2 














e-i p p p 


P 


«x^, 


P P 














p -i/> •-•W-CO- 


• 


Eh VV 














P «r H ^ «c 


m 


P 


<r:^ 














P Eh CL, 6-1 £h 


c!j 


P 


£-< P Eh 


• * 


• » 










P 2 p P 


P 


P 


P2P 


— • 












1 P p W p 


P 




P P P 


to 












6-< P P w 


+ 


P Eh P 


tH 










4t 


^ • • • • 


O 




• • • 


w 


-iO- 




■w* 




-> 


P • • • • 


Eh 


e? 


» » * 


o 


C!3 




•Ji* 




4c 


P O ci? P c3 w 


p 


OOP 


CO 


CO 




*/% 




4$- 


2 P P 2 2 


P 


p 


P p 2 




2 — 




*K- 




4C- 


Z ffi p p p 


«x. 


2 


P P P 


• 






•K* 




■J?- 


i-» 2 2 E-* 2 


P 


2 


2 2 Eh 




— X ••- 




<- 




-X- 


CO z CO • • 


• 


HH 


2 CO • 




X ^ 




-5^- 




4C- 


»— 1 o • 


•> 


CO 


o • 


o 


C: X 




•K* 




45^ 


45* CO O Ci) o 






CO C-) o 


2 


2 • X 




rc 




4C- 


• • Qd P 


P 


45* 


• • 2 




*-H •• pq 








4i- 


^ • *p p 


P 




• •* P 


a: • 


te 9 h • 




•ic- 




45* 


c >-• xwh;>“hA* 


X 


X >hha 


E-i^ 


E-4 




<- 




45* 


6-, <r <i; ^ <1. 


«c 


«r, 


<C 


oo 


CO 




•VC 




4c 


p 6- Eh Eh 6-^ 


Eh 


Eh 


Eh Eh 6- 


w cti 


•w X cq 




•K* 




45- 


p p p p p 


P 


P 


P P P 


6h O 


Eh X O 

OTj O" 








45* 

>5, 


p p p p p 
r~\ 


P 


P 


P P P 




^ P-J M- 








•5c 


— ^ P P P P 


P 


J— H 


P P ^ 


i-H p:^ pL, 


*-H • a: 




•K- 




45* 








• • • 


p- W »-q h:* 1^ 


pc -- W 








45* 


II > ''w- 




ii 


W 'w-' 


W cc pt- oc 


P -wn- 


p 






45* 










W o o o 


-^A »-:? P< 


o 






45* 


c3 P P P >* 


P 


O 


P P P 


£h CO i/>'</> W 


6h CO 


-ify 


<- 




4t 


p :d :=> ♦-» 


P 


2 


POP 


pcj »-i n 


pp 5: »-H 


n 






45* 


P 2 2 CO P 


2 


P 


2 2 <r 


O P^ 2 2 2 


O f*. P5 


2 


45- 


CO 


45- 


w P P P P 


P 


-lA 


P P P 


P*J Pr' W 




W 


•i'f 


p; 


45* 






< 




►-a CO CO cn 


►J »-:» »-4 


CO 


•Jc 


o 




6h p p p p 


P 


EH 


P P P 


l-J 1-^ 


►J »-:» 1-^ 




4C* 


pcJ 


-X- 


P P p P P 


P 


P 


P P P 


1-^ »-3 1-5 


«c *< 


P5 


4c 


P 


4c 


p <i: <t: <j^ 


<r 


P 


< «x^ <r 


O O h 3 1-^ l-J 


o o o 




*rC 


P 


45* 


p o o o o 


o 


P 


O CJ) o 


<C -or, «< 




<c 


4# 

/4 




45- 










O o o 




O 


4i- 




45* 


45* 




45* 










\4t 


4^ ^ 


45* 






\ 





C\i 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CO 




m 




^- 


cc 


Gi 


s: 


P 


CM 


CO 




Ml 


CO 


l> 


00 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CM 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


C^J 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 








cc; 




\ 




2: 








\ 




















cj 




pc; 




Pi 




PQ 




pp 




2: 


• m 


z: 




CO 




CO ••• 




0 


EH 


0 ^ 




0 


0 


0 EH 






pp 


0 






H/> 


^ PP 






0 


• •* H/V 






Pi 


Eh^ ^ C 3 




0 


EH pp 


0 6 h £h 




0 


0 w 


no 0 pi 




-uv 




W pp pp HTV 




X 


W ptl 


>-4 -o> W hP 




••* 


0 PP 


hP 0 0 W 






« ••- •• • 


z: •- 2: PP 




rP EH E-» 




n Cd — pH • 




0 0 


V> Pm P-* Eh 0 


-W PhVV •“ 




1 0 


hP 2: z: 0 w 


+ hP z: i-p 0 




V> 1 /V 


Cxi w PM C 5 PP 


PM PM P«q 




0 X >H 


PP £h Eh VV W 


0 PP Eh PH CO 




hP hJ 


• • • c 3 0 


pH*«*vV 




PP pd 


•‘(Xi Eh 


pq - 0 




2: P) PP 


c!j 0 P- P^l H/V 


z: 0 0 PM EH 




Z • • 


pc: pd ziha pp 


2: Cd PM 2 T -o> 




H- 4 ^ • 


PQ pp PM ,-JI < 


pp pq PM Pd 




CO PP PP 


2: 2: Eh fxl p- 


CO 2: Z. EH PP 




EH Eh 


z: CO • PP • 


3: CO * * 




■ 5 t CO CO 


♦— 4 0 • • • 


hh 0 • • 






CO 0 6 -< - Eh 


CO 0 Eh EH 




EH ►<; >-. 


• -COO 


Eh • • 0 0 




0 > >• 


- •‘CP Z PP 


0 •‘PP PP 




0 * • 


6 -^ Eh H/> Q- -cO* 


n Eh EhhA hA^ 




•V> • • 


0 0 e z: e 


hA 0 0 0 0 




>H «C ^ 


CP pp pr, •Pi 


X PP PP 2 3: 




hP EH EH 


V>H/> pp •* pp 


hPhAha pd Pi 




W CO CO 


XX-«/>PMh/> 


C*i X Xv> ha 




PP P«^ w 


hp iJ z: 1-5 


PH hP n n i-M 




^ >M 


pc. PM w Pm 


PM PM P*^ PM 




II >► >• 


pp PP P; Eh PP 


)i PH n PP pH 




• • 


• • • • • 


• • • • 


CV* 


EH w- — 





g-i — - 




0 




0 




PP PQ 


►-q hP pr. hP 


P i-p n n hp 




H/V 5 D 0 


pii SD CD »-H :d 


hA p p Pi P 


w 


0 CO CO 


z: 2: CO pp 2: 


e 2: z- 2: 




pc: P*^ . Pm 


Cm Ph pM Pm 


2: Pm pt 4 pu Pm 




pn 




P- 


P- 


H/V 1 ^ hP 


,J iJ hJ hP hP 


W p p p p 




I-P hP i-P 


i-P iJ »-P hP 


p p p p p 




P«a -11 < 


< 1 . <1 «a. . 


P -d -t. «tj 




pp 0 0 


0 0 CD 0 0 


p 0 0 0 0 


pci 


*«• 




<• 


W 






\ 










H - 1 








PM 








5: 








0 








0 










CVi c\2 


CVJ CVi CJ CVI CVi 


0 ^ CM CVi CO 


0 








CD 

1 








1 

z; 


C 7 > <Si 


irH CVi to m 


0 cc cn 


\ 


CO 


t}^ ^ 


kH ^ 




CO CO 


CO CO CO CO CO 


CO CO CO CO 


PM 









TOT 



• 








• •k 


















V 






• «k ••k 








^ — . 




m 








CD 














cs 






















to — 




tH 








CD 


Eh 




f-H 






a> Eh 






If) 






w 








HA 








ha O 




hA 






hA O 






HA 


O C!) 




c!> 




C5 c> 




C5 








O 


P 




O 






C!) P 






o 


in z 




to 




to P 




to 








CO HA 




to 






tOH/^ 






in 


z: tti 




Z •- 




ZP 




z 


• •> 






Z 


O 




z 






zo 






z 


• -»/> 








• w 




• 








• 


P 




« 


V 




• z 






« 


' — ^ «c 








— ^ <d 




'^- 


X 









P 




' — 












— ' 


E-* 




X-— 




e-^ 






X 








ha 










•Ur 








C3 h:^ 




O P X 




O P 




ci? 


P X 




o 


p 




C' 


P X 




c!) P 






o ••* 


z w 




Z • X 




Z P 




z 


« 


X 




z 


W 




z 


• X 




z W 






z o 


1-^ pi 




l-H » p 




•-H p 




»-H 




p 




»-H 


P 




»-H 


•p 




HH pi 






»-H Z 


a: • 




cci e-. • 




p • 




p 


Eh 


« 




p 


• 




p 


Eh • 




P • 






PhA 






e-* • ^ 








e-t 


• 


' — ' 




Eh 






Eh 






Eh'w^ 






Eh to 


in 




in • 




in 




to 


* 






to 






CO 






to 






to W 


w cci 








ha « 




H/> 


X P 




hA 


P 




HA 


X P 




hA pH 






HAX 


e-« o 


« •> 


E-^ X o 




e-< o 


• m. 


E-* 


X o 




Eh 


o 


• 


Eh 


X o 




Eh O 


• •> 




EhhA 


z: cc P»4 


P P 




P 


P z 


p 


p 


• •» 


Z 


p 


w 


Z 


p p 




Z p W 


W 


z; 


»-• ft: 1-4 


»-H . Pi P l-H p 


p »— 


• 


p 


W *-i 


p 


p 


KH 


• P w 


i-H p p 


P 


»-H O 


oc: fk: pc: 


Pd^ P h4 p P 


cc 


p: 


' — 


WPP 


w 


et. 


P 


W P 


P W 


P 


P 


P w 


Ph hA 


c; 


p hA P 


O P 




-a> P P 


ha 


o 


P 


HA P 


PhA 


O 


o 


p. m 


w 


w 


-U> i-a P O-lA P 


-0>H/h 


p 


P OhA 


P 


HA 


hA 


P P o 


hA P- 


hAHA 


hA O 


e-» cQ 


Pi 


E-* CO W 


£-• CO 


P 


£-• 


p 


CO 


HA 


Eh 


to 


P 


Eh 


rz) CO 


HA 


Eh CO 


p 


p 


E^ 


P=i 


z 


P Z 


p 


P HH 


z 


P 


z 


»-H 


p 


P 


»-H 


Z 


P 


Z »-H 


p 


P »-H 


Z 


z 


p. n 


o p 


p 


O P p 


z 


o p 


w 


O 




P 


Z 


O 


P 


w 


O 


W P 


z 


o p 


W 


w 


c. 




00 




w 




to 








w 






to 






p 




CO 


CO 


X 


1-4 P 




h4 P P 


to 


PP 




P 


p 


P 


CO 


P 


P 




P 


p p 


to 


p p 






p <«; 


1-4 1-4 


h4 


»-4 P 1-4 




P P 


p 


P 


p 


P 




P 


P 


P 


P 


p p 




p p 


p 


p 


p P 


<<r 




^ <£J 


p 




p 


«r. 


<€t; 


•< 


p 


«c 


•«<; 


P 






p 


<t. 


p 


p 


«-< p 


O O 


<Jh 


o o o 


p 


O O 




O 


O 


O 


p 


o 


o 




CJ 


o o 


p 


CJ) O 




<< 


O 




O 




< 




o 














o 










o 


o 










o 












o 












o 











CVi 


cv 


CV2 


CVi 


CVi 


CV2 


04 


04 


CVi 


04 


04 


C4 


04 


04 


04 


04 


04 


C4 


OJ 


C4 


04 


04 


04 


04 


04 


CV) 


04 


cs 


p 


CV2 


CO 




If) 


CD 


0- 


CO 


0) 




P 


04 


CO 




m 


CD 


C^ 


00 


O) 




P 


04 


CO 




m 


CD 


ifj 


If) 


\il 


m 


If) 


If) 


If) 


If) 


If) 


in 




to 




CD 


CD 


CD 


CD 


CD 


CD 


CD 


0- 




0- 


!>- 


O- 


c- 


£> 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


CO 


C^) 


co 


CO 


CO 


CO 


CO 


CO 


CO 



128 



377 2 I = I + i; 

378 2 ?ND; 

379 1 END MAIN^PROG; 
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NUMOUT: 

PROCEDURE USED TO CONVERT A NUMBER FROM (INTERNAL) BINARY 
FORM TO AN ASCII STRING SUITABLE FOR CRT. IT IS WRITTEN 
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riSPLAY^CRT ; 

PROCEDUPP USED TO DISPLAY ANY POUR “BYTE VARIABLE IN THE 
EORM OE 5 INTEGER AND 2 DEGIMAL PART OE IT. 
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IF CHAR >= 80H THEN CHAR = CHAR XOR 80H 
RETURN CH*^R; 

END crt$read; 



CHFCK^YTS^NO : 

PROCEDURE USED TO CHECK FOR A VALID YES/NO INPUT FROM THE CRT 
US.*GF ; TYPED PROCEDURE, ft VftLUE OF 1 WILL PE RETURND IE THE 
ANSWER IS YES OTHERWISE THE VALUE RETURNED IS 0. 
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120 2 END DISP^ERROR; 

121 1 END display; 



APPENDIX F. FLOATING-POINT PROGRAM LISTING 
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DO WHILE (DUFFER (POINTER) O EOL ) OP ( BUFFER ( POI NTER 
CALL CRT$WRITE (BUFFER ( POINTER )) ; 

POINTER = POINTER + IJ 

fnd; 

END crt$print$string; 
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0UTPUT(MEM$L0W) = LOW(MEBAS); 
OUTPUT(MEM$HIGH) = HIGH(MEBAS); 

return; 

END INITiFP; 
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0P2(I) = PES*T ABT.Ed ); 
OPKI ) = RES$TA.BLE(I + 

end; 

fetupn; 

END VAL$RESULT$i; 



VAL$RESULT$2: 

* PROCEDURE USED TO GET THE RESULT FROM FIXED POINT 
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CO^^PAHE: 

* PROCEDURE USED TO CHECK EOR OUTPUT CONDITIONS FROM E.P. BOARD 
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C/iLL CRT$PRINT^.STRlNGf .MSGl 
CALL CRT^PRINT$STRING( .HSGP 
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CALL ABJUST^OP (<^,P); 

IF (RES$TABLE(3) >= 80H ) AND 
(PES*TABLE(7) >= 80H ) 

THEN FLAG = 0FFH; 
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END 



end; 

CONVERT ANGLE IN RADIANS TO ANGLE IN SEMICIRCLE UNITS 
CALL FDIV( .TEMP 0 , .PUnO^T , .TEMP); 

GET THE SQUARE OE THE ANGLE. */ 
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GIVE SIGNS TO COSINE AND SINE VALUES ACCORDING TO QUADRANT 
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ANGLE(I) = PU3$0VEP.2(I ); 

end; 

return; 



I = 0 TO 3; 

ANGLE(I) = PI$0VER2(I ); 

end; 
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CALL FADD(.TEMP, .C0NST$1, .TEMP); 
CALL FMUL(.Z, .TFMP, .TEMP); 

CALL FADD(.TFMP, .PI$0VTi'R4, .ANGLE); 
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